Dirac returns: non-Abelian statistics of vortices with Dirac fermions 
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Topological superconductors classified as type D admit zero-energy Majorana fermions inside 
vortex cores, and consequently the exchange statistics of vortices becomes non-Abelian, giving a 
promising example of non-Abelian anyons. On the other hand, types C and Dili admit zero-energy 
Dirac fermions inside vortex cores. It has been long believed that an essential condition for the 
realization of non-Abelian statistics is non-locality of Dirac fermions made of two Majorana fermions 
trapped inside two well-separated vortices as in the case of type D. Contrary to this conventional 
wisdom, however, we show that vortices with local Dirac fermions also obey non-Abelian statistics. 
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One of remarkable achievements in condensed matter 
physics in these years is the classification of topologi- 
cal insulators and superconductors [l|, 0). In particu- 
lar, among them, the type D superconductors exhibit 
a unique property that zero-energy Majorana fermions 

H| are trapped inside a vortex core @ or appear at the 
edge of superconductors. Examples of materials for type 
D are given by chiral p-wave superconductors/superfiuids 
such as Sr 2 Ru04 and the A-phase of 3 He. One of im- 
portant aspects of the zero-energy Majorana fermions is 
that exchange statistics of vortices is non-Abelian 0, 0] 
so that a set of them can be regarded as non-Abelian 
anyons. They are expected to be quite useful for a fault- 
tolerant quantum computing 0] . The idea for topological 
quantum computation utilizing non-Abelian anyons on a 
topological state of matter has been developed [§] . This 
is why it is one of hot topics to experimentally realize 
Majorana fermions. 

In the first theoretical discovery of non-Abelian statis- 
tics of vortices with Majorana fermions (we call them 
the Majorana vortices) [7|, the key ingredient was an 
unusual way to define Dirac fermions that are neces- 
sary for the construction of a Hilbert space of many 
vortices. With a single Majorana fermion at each vor- 
tex core, one has to compose a Dirac fermion from two 
Majorana fermions located at spatially well-separated 
vortices. Since then, it has been long believed that 
the essential property that leads to non-Abelian statis- 
tics is the non-locality of Dirac fermions. This is one 
of the reasons why other types of topological insula- 
tors/superconductors have not attracted more attentions 
than the type D insulators/superconductors as candi- 
dates of non-Abelian anyons. 

On the other hand, the type C and Dili insula- 
tors/superconductors admit Dirac fermions inside vor- 
tex cores [Io| . We call such vortices the Dirac vortices. 
There are several systems allowing Dirac fermions such 
as cores of integer (singular) vortices in 3 He A-phase 



and dislocation lines in topological insulators [12| . Dirac 
fermions already exist locally in these cases, and there- 
fore they have been believed for a long time to admit 
no non-Abelian statistics. Contrary to conventional wis- 
dom, however, we show in the present paper that they 
actually do give non-Abelian statistics, and consequently 
it implies that Majorana fermions are not necessary any- 
more for constructing non-Abelian statistics. 

We treat a set of n vortices, each of which allows for a 
single zero-energy Dirac fermion at its core, and consider 
an exchange of two vortices in such a system. An opera- 
tion T k is defined as the anti-clockwise exchange of the fc- 
th and (fc+l)-th vortices. Then, {Tk\k = 1, • • • , n— 1} are 
in general the generators of the braid group. They satisfy 
the braid relations (i) T k T e T k = T e T k T e for |fc - 1\ = 1 
and (ii) T k T t = T e T k for \k - £\ > 1. We assume that 
exchange of two vortices is performed as an adiabatic 
process, and that quantum states of vortices are essen- 
tially determined by the zero-energy modes. Let us define 
an (annihilation) operator tp k for the zero-mode Dirac 
fermion at the fc-th vortex. It satisfies the following al- 
gebra: {V>fc,V>j} = 5u, {4>kAl} = = 0. Then, 
the exchange of two vortices is expressed in terms of the 
change of operators ip k under the operation T k . Note that 
the wave function of the zero modes obtains an additional 
phase, a minus sign, when the zero-mode fermion turns 
around a vortex, implying that the zero-mode wave func- 
tion is a double- valued function. In order to regard this 
wave function as a single-valued function, we introduce 
a cut so that the zero-mode wave function changes its 
sign when it crosses the cut. Under the operation T kl we 
find that the fc-th vortex crosses the cut of the [k + l)-th 
vortex. The operation T k acts on the Dirac zero-mode 
operators ip k and i/'fc+i for the fc-th and (k + l)-th vor- 
tices, respectively, as 



Ti 
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with keeping the others ipe (£ ^ k,k + 1) unchanged. 1 
We confirm that the above transformation satisfies the 
braid relations (i) and (ii). We should emphasize that, 
because of the minus sign in Eq. ([1]), we obtain the re- 
lation Ti = 1, i.e., four successive exchanges are equiv- 
alent to the identity. 2 This may be contrasted with the 
exchange of two conventional particles, in which two suc- 
cessive exchanges are equivalent to the identity. 3 This 
fact by itself suggests that the operation Tk induces a 
non-trivial statistics, as discussed below. 

Now, let us find the representation of the operation Tk 
following the procedure discussed in Rcf. . First of all, 
we note that Tk defined in Eq. ([T]) can be represented 
with respect to the Dirac fermion operator tpi. The cor- 
responding operator Tk reads 

Tk = 1 + ipk+i^l + ^fc+i^fc 

-^fe+i^fe+i ~ W^fc + ^l+i^fc+i^IVV ( 2 ) 

We confirm that Tkinf^ 1 {I = 1, . . . ,n) reproduces the 
transformation ([!]). We also note that fk satisfies the 
braid relations: (i) f/cffffc = t^t^t^ for |fc — £\ — 1 and 
(ii) ffcff = T£Tk for \k — l\ > 1. Moreover, we have 
fi = 1. Therefore, fj. gives a representation of the braid 
group with a condition = 1. In general, f^'s are non- 
Abelian; TkTk+i ^ Tk+\Tk, which will be explicitly seen 
in the matrix representation of T£ in the Hilbert space. 5 

Next, we note that, with the operators ipi {I = 
1, • • • , n) for the zero- mode Dirac fermions, it is straight- 
forward to construct the Hilbert space of the multiple 
vortex systems. The Fock vacuum |0) is defined by 
V-v|0) = for all £, and one can construct the Hilbert 
space by multiplying ^t's on the vacuum. Since the 
fermion number / is conserved under the transformation 
@ (the fermion number operator / — ^"=1 com ~ 
mutes with fj, for any fc), the total Hilbert space M^ n ' 
for n vortices is decomposed into a direct sum of the 
sectors H^™'^ with the fermion number / (0 < / < n); 
ffl(™) = ®^ =0 By successive multiplications of ip\ 

on |0), we obtain states with the fermion number /. We 



1 The transformation JTJ is consistent with the Bogoliubov-de 
Gennes equation which is a microscopic equation for the fermions 
in the presence of vortices. 

2 Two successive exchanges yield T 2 = — 1. 

3 If there is no minus sign in Eq. (£TJ , we obtain = 1 which im- 
plies the Bose-Einstcin/Fcrmi-Dirac statistics or para-statistics. 
Thus, the existence of the minus sign in Eq. JTJ is essential. 

4 Recall that, in the case of the Majorana vortices, four-times 
multiplication of an exchange operator is projectively equivalent 
to the identity [3]. 

5 In contrast, f^'s are Abelian when fj.'s are just complex numbers 
e' t with real 9^'s (9 = for bosons, 9 = n for fermions, and 
9 ^ 0, 7r for anyons). 



use the following notation: 

|0...01...1...10...0)=4V--4V--$ / l°>> ( 3 ) 

1 h k l f n 

where l\ < ■ ■ ■ < t+ < ■ ■ ■ < If, 1(0) in the i-th for the 
(non-)existence of the t-th. Dirac fermion in the left hand 
side. The /-fermion sector M^ n '^ is spanned by these 
states and has dimension n Cf = yr^zyyr- 

The matrix representation of the operator Tk is also de- 
composed into n Cf x n Cf blocks Tj£ n , corresponding to 
the /-fermion sectors M.' n, f\ For example, for the basis 
with and/or 1 for fc-th and (A; + l)-th Dirac fermions, 

the elements of the matrix t^™'^ are given by 

...oo... h ...on... = l, 
...io... h ...oi ... = -l, 
...oi ... h ... io... = 1, 

(...ll...|f fc |...ll...) = 1, (4) 

and the others are 0. Notice that the matrix representa- 
tion is allowed for both even and odd number of vortices. 6 
Then, the total matrix is given as = (B/=o r i™ ^- 
We will see below that this basic structure essentially de- 
termines the non- Abelian statistics of the representation. 

Let us see concretely how to construct the Hilbert 
space M(") and how the non- Abelian structure emerges 
in the matrix expression of t^™ (0 < / < n) by using 
simpler cases with n = 2,3,4 vortices. In the case of 
n = 2, we have only one operation T\. The Hilbert space 
Hi' 2 ) can be decomposed into a direct sum of three sectors 
H( 2 <°) = {|00)}, rf 2 ' 1 ) = {|10), |01)} and rf 2 ' 2 ) = {|11}}. 
Then, we obtain the matrix representation of f\ as (see 
Eq. Q) 

# 0) = 1, ^-(^o 1 )- 4^=1. (5) 

Here, t[ ' and t[ ' for the cmpty/fully-occupied sec- 
tors H( 2 <°) and H( 2 ' 2 ) are identities, while rf' 1 ' for H^ 1 ' 
contains off-diagonal components. Even though the ma- 

(2 1).. 

trix t± ' is non-diagonal, we have only one such matrix, 
and thus can simply diagonalize it to find an Abelian 

(2 1) 

statistics. Indeed, ' can be diagonalizcd by a uni- 
tary transformation, and we find that it has eigenvalues 
±i. Therefore the exchange statistics is an Abelian anyon 
statistics in the sector rf 2 ' 1 ). 



It should be noticed that it makes sense to discuss the exchange 
of occupied "1" and empty "0" states which corresponds to the 
matrix elements, such as (. . . 10 . . . |ffe| . . .01 . . . )■ This is be- 
cause we are considering the exchange of two vortices, and it is 
possible to have the situation where the vortex does not contain 
a Dirac zero-mode fermion in it. 
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Truly non-Abclian statistics emerges when the number 
of Dirac vortices is larger than or equal to 3. In the case 
of n = 3, the Hilbert space H( 3 ) can be given as a direct 
sum of H( 3 >°) = { 000)}, H^' 1 ) = {|100>, |010), |001», 
H( 3 > 2 ) = {|110>,|011>,|101>} and rf 3 ' 3 ) = {|111>}. We 
thus have the matrix representation of fk as 



.(3,0) 



.(3,1) 



.(3,2) _ 



.(3,3) 



.(3,0) 




.(3,1) _ 



.(3,2) _ 




(0) 



The empty/fully-occupied sectors H^ 3 ' ^ and H^ 3 ' 3 ^ are 
again trivial. For the sectors IH^ 3,1 ) and H^ 3,2 ) the ma- 
trices show the basic structure ((3]). This time, it is not 
possible to simultaneously diagonalize all the matrices. 



,/) 



and r 9 



,/') 



are non-commutative; 



Indeed, we find t{ 

r (»./) r K/) ^ T ^" f) r[ nJ) with n = 3 and / = 1,2. 
Therefore the exchange statistics is non-Abelian in the 
sectors H^ 3 ' 1 ) and H^ 2 ). 



We again note that, with the Dirac fcrmions, the 
Hilbert space is constructed for both even and odd num- 
bers of vortices. This is in contrast with the vortices with 
Majorana fermions 0, in which case the Hilbert space 
can be constructed from only even numbers of vortices 
because a single Dirac fermion is constructed from two 
Majorana fermions belonging to different vortices. 

Finally, let us also show the results for the case of 
n = 4. We will see the basic structure (0| and thus 
non-Abelian statistics again when the fermion number 
/ = 1,2,3. The Hilbert space rf 4 ) can be given 
as a direct sum of five sectors: H^ 4 ' ) = {|0000)}, 
H' 4 ' 1 ) = {|1000>, |0100), |0010), |0001)}, H( 4 ' 2 ) = {|1100>, 
|1010), 1001), |0110), |0101), |0011)}, M^ 3 ) = {|1110), 
jllOl), jlOll), |0111)}, and H< 4 - 4 ) = {|1111)}. The ma- 
trix representation of fk is: 



.(4,0) 



.(4,1) 



.(4,1) 



.(4,0) 
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in the sectors H^ 4 ' ) and rf 4 ' 
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in the sector H' 4 ' 2 \ and lastly 



.(4,3) 
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in the sectors IH^ 4 ' 3 ) and H( 4 ' 4 ). As in the previous cases, 
the empty /fully-occupied sectors H^ 4,0 ^ and H^ 4,4 ' are 
trivial. For the sectors H^ 4 ' 1 ) H^ 4 ' 2 ) and U^ 3 \ we find 



JnJ) (n,f) 

T l ! T 2 

JnJ) (n,f) 
T 2 T l > 



and 



>,/) 



JnJ) (nj) 
T 1 T 3 



non-commutative; t\ 



(n,f)(nj) 



+ r. 



(n,f)(nj) 



with n = 4 and 



/ = 1,2,3. Then the exchange statistics is again non- 
Abelian in the sectors H^ 1 ), H< 4 ' 2 ) and rf 4 < 3 ). In this 
way, we can construct the matrix representation of the 
exchange operators ffc for arbitrary number of Dirac vor- 
tices. 

Our procedure to find the matrix representation for 
the exchange of two Dirac vortices is very similar to the 
one adopted by Ivanov 0] for the exchange of two Majo- 
rana vortices. However, let us clarify that these two are 
quite different in several aspects and emphasize that the 
complexity appearing in the manipulation of the Majo- 
rana vortices is not an essential requisite to realize the 
non-Abelian statistics. We here point out two aspects: 
locality of the Dirac fermions and conservation of the 
Dirac fermion number under the exchange of two vor- 
tices. While these two hold for the Dirac vortices as 
explicitly demonstrated in the present paper, they are 
violated in the Majorana vortices. 
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Concerning the locality of the Dirac fcrmions, we re- 
call that the Dirac fermions employed to construct the 
Hilbert space for the Majorana vortices are nonlocal. 
More precisely, the Dirac fermion operator ip^f = \ {^k + 
il2k+i) is made of two Majorana operators 72/c and 72fc+i 
in spatially separated 2fc-th and (2k + l)-th vortices [7J. 
Namely, 2n Majorana vortices generate n nonlocal Dirac 
fcrmions. In contrast, since we started with the Dirac 
vortices themselves, they are by definition local, i.e., a 
Dirac fermion is located at the core of a single vortex. 

In the Majorana vortices, non-Abelian statistics ap- 
pears when one exchanges two 'unpaired' Majorana vor- 
tices, each of which is used to define different Dirac 
fcrmions. Then, the exchange of two unpaired Majorana 
vortices induces non-conservation of the Dirac fermion 
number Q. For example, exchange in the zero Dirac 
fermion sector generates a state with two Dirac fermions. 
Mathematically, one finds that the fermion number op- 
erator / M = J2k=i V'fe^V'fe 1 which counts the number of 
Dirac fermions does not commute with the exchange op- 
erator t%1 +1 = -4= (1 + 72/c+272fe+i) corresponding to the 
exchange of the unpaired Majorana vortices. In contrast, 
as mentioned before, the fermion number in the Dirac 
vortices is conserved under the exchange of two vortices. 

We find interestingly enough that the Dirac operator 
in Eq. J2]) can be decomposed into two Majorana opera- 
tors fjl and fjl = f)l) as tp k = \ (f\\ + ifj 2 k ) , up to an 
overall phase. (Notice that this expression is consistent 
with the exchange rule ([T]), because each Majorana oper- 
ator rf k changes the sign like Eq. ([U under the exchange 
of vortices.) Then, the operator ([2]) can be rewritten as 
f k = Ili=i,2 73 ( l + ijl+iVl) ■ Tnis suggests that if wc 
start with two Majorana fermions localized in a single 
vortex, then the Hilbert space of the Majorana vortices 
can be spanned by local Dirac fermions, and the fermion 
number is conserved under the exchange of vortices. 

However, if we have three, or in general, odd num- 
bers of Majorana fermions in a single vortex, we can- 
not construct Dirac fermions only by using the Ma- 
jorana fcrmions in a single vortex, and thus necessar- 
ily encounter non-locality of Dirac fcrmions and non- 
conservation of fermion number. In fact, this was ex- 
plicitly demonstrated when the number of Majorana 
fermions is three and they have SO(3) symmetry [l3j |. 
The resulting non-Abelian statistics was found to be 
given by a tensor product of the matrices in a single Ma- 
jorana case and a Coxeter group. An example of physical 
systems to realize it is a color superconductor [14j. 

Hence, we observe that non-Abelian statistics can be 
realized by using either even or odd number of Majo- 
rana fermions. When the fermion number is even (odd), 
one uses local (nonlocal) Dirac fermions, and the fermion 
number operator is (is not) conserved under the ex- 
change of two vortices. Therefore, we emphasize that 
both the non-locality of the Dirac fermion and the non- 
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FIG. 1: The exchanges of (a) two sets of two Majorana 
vortices and (b) two Dirac vortices, which yield the statistics 
given by Eq. (JTUJ and Eq. (TTJ, respectively. They are not 
equivalent; the latter is non-trivial statistics studied in this 
paper, while the former is trivial statistics. 



conservation of the fermion number are not essential req- 
uisite for the non-Abelian statistics, but the minus sign in 
the exchange of vortices as in Eq. ([1} is rather important. 

One interesting question is if two Majorana vortices 
behave similarly to one Dirac vortex, and moreover if 
exchange of two sets of two Majorana vortices show 
non-Abelian statistics. However, the exchange of two 
Dirac vortices is not equivalent to exchange of two sets 
of two Majorana vortices, see Fig. [1] Consider two 
Dirac fermion operators in a system of (even number 
of) the Majorana vortices: tp^f = 5(72*; + «72fc+i) and 

4>Wi = ^(72fc+2 + *72fc+3), where the first one is made 
of the 2fc-th and (2k + l)-th Majorana vortices, and the 
second one the (2k + 2)-th and (2k + 3)-th Majorana 
vortices. Using the basic transformation rule for the Ma- 
jorana fermions, one can easily find the exchange trans- 
formation of these two Dirac fcrmions: 



Yk+1 



4 M 



(10) 



Thus, there is no minus sign, unlike Eq. ([TJ. Sec Fig. [1} 
(a). The matrices for exchange operators fk in this case 
can be obtained by ignoring the minus signs in our ma- 
trices. It satisfies fj? = 1 giving at most parafermion 
statistics, see footnote [3] 

Finally we make a comment on supersymmetric gauge 
theories as a possibility of Dirac vortices. Bogomol'nyi- 
Prasad-Sommerfield vortices spontaneously break and 
preserve a half of supersymmetry in supersymmet- 
ric gauge theories. Consequently Nambu-Goldstone 
fermions (Goldstino) for spontaneously broken supersym- 
metry is trapped in the core of the vortex [l5|]. When 
the theory posseses M = 2 supersymmetry (eight super- 
charges) M = 1 supersymmetry (four supercharges) is 
preserved. Therefore the fermion trapped in the vortex 
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is a Dirac fermion because it belongs to a chiral multiplct. 
Therefore vortices in supcrsymmctric theories may give 
non-Abelian statistics studied in this paper. It remains 
as a future problem to study this case. 

In summary, we have explicitly constructed the ex- 
change statistics of vortices with Dirac fermions, which 
is relevant to the type C and Dili topological insula- 
tors/superconductors and some other systems, and have 
found that it is non-Abelian, contrary to conventional 
wisdom. It is expected that the existence of local Dirac 
fermions will provide us an interesting approach to ob- 
serve the non-Abelian statistics in laboratories, helping 
us to utilize quantum computations. 
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